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ABSTRACT 



The next-to-leading order QCD corrections to the production of gluino pairs at the Tevatron 
are presented in this paper. Similar to the production of squark-antisquark pairs, the depen- 
dence of the cross section on the renormalization/factorization scale is reduced considerably 
by including the higher-order corrections. The cross section increases with respect to the 
lowest-order calculation which, in previous experimental analyses, had been evaluated at the 
scale of the invariant energy of the partonic subprocesses. 
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1 Physical Set-Up 



The novel colored particles in supersymmetric theories, squarks and gluinos, can be searched 
for most efficiently at high-energy hadron colliders, i. e. the Tevatron pp collider and the LHC 
pp collider in the future. The most stringent lower bounds on squark and gluino masses have 
been set by the Tevatron experiments CDF and DO. At the 90% CL, the lower limit of the 
gluino mass was found to be 157 GeV, independent of the squark mass [[[], @- 

The cross sections for the production of squarks and gluinos in hadron collisions were 
predicted at the Born level already quite some time ago Only recently the predictions 

have been improved by the next-to-leading order QCD contributions for squark-antisquark 
pair production In a first step, the gluon radiative corrections to gluino pair production 
have been calculated in Ref . || ; these corrections are closely related to the gluon corrections 
for quark-antiquark production, requiring just the appropriate change of the SU (3) Casimir 
invariants. In this paper we shall present a complete next-to-leading 0(a s ) analysis for 
gluino pair production, 

pp^gg + X (1) 

including the virtual effects of gluinos and squarks besides gluons and quarks. 

The technical set-up of the calculation is analogous to the squark analysis. For the sake 
of simplicity all squark states are taken mass degenerate. The only free parameters are 
therefore the masses of squarks and gluinos, rriq and rrig, respectively. (The top mass is fixed 
to 176 GeV g [nj.) 

The calculation has been performed in the Feynman gauge and the singularities have been 
isolated by means of dimensional regularization. The masses are renormalized in the on- 
shell scheme. The massive particles are decoupled smoothly for momenta smaller than their 
masses in the modified MS scheme [JTTJ] . The soft and the hard gluon radiation are separated 



by a small invariant mass cut-off parameter A for the gluino-gluon and gluino-light quark 
final states, regularizing the infrared divergences after the virtual corrections are included. 
If soft and hard gluon emissions are added up, the A dependence disappears from the total 
cross section for A — > 0. The remaining collinear mass singularities are absorbed in the 



renormalization of the parton densities carried out in the MS scheme [|12] . We have adopted 



the GRV |13| and CTEQ(2pM) |14| parametrizations of the parton densities to compare 
the NLO predictions with the LO approximations; the MRS(H) |15j parametrizations have 
been included finally to assess the uncertainties of the predictions associated with different 
parametrizations of the parton densities. 

If the gluinos are lighter than the squarks, on-shell squarks can decay into quarks plus 
gluinos, q — > qg. We will assume in the present analysis that these events are attributed to 
squark final states; technically, we subtract off the squark-gluino production cross section 
Q9 ~^ Q9 after including the non-zero squark width in the propagator. Analogously, the 
decay of gluinos to squarks is disregarded in the wedge m- g > nig. Spurious singularities in 
higher orders at the decay threshold for gluino decays to stop and top, can be regularized 
by taking the non-zero widths of the particles into account and by performing the mass 
renormalization of the gluino at the complex pole of the propagator. 

The dominant contributions to the production cross sections of gluino pairs in pp collisions 



2 



are due to the gg — > gg and the qq — > gg subprocesses. In contrast to squark-pair production, 
the gg initial state is dominant if both m g and nig are less than about 200 GeV. The qg 
initial states give rise to gg pairs only at next-to-leading order while the contribution of 
quark-quark pair initial states to gluino pairs are of yet higher order. 



1.1 Gluon-gluon initial states 

To lowest order the diagrams contributing to the subprocess gg — > gg are analogous to the 
diagrams for quark-pair production, Fig. la. Typical standard QCD and supersymmetric 
vertex corrections are displayed in Fig. lc. Gluon radiation processes, Fig. Id, are added 
incoherently. 

Evaluating these diagrams analytically we obtain the double differential cross sections 
d&ij j dtdu at the parton level; i,j are the parton indices g,q,q and t,u the Mandelstam 
momentum transfer variables. The total cross sections may be expressed in terms of scaling 
functions fy, 



al{Q 



a t3 = ^^{flf + Ana s {Q 



m g 



+ =2) } (2) 



Q 




The functions depend on the invariant parton energy and the gluino/squark masses. 
a s is the QCD coupling constant. For the sake of simplicity, the renormalization and factor- 
ization scales are identified, Hr — ^f — Q- For the non-zero scaling functions to lowest 
order compact expressions can be derived, 



f(0) = ^ 
■I 99 5 
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where (3 — (1 — 4m~/s) 1 / 2 denotes the gluino velocity in the parton c. m. system. The scaling 
functions f-p and f-p, which describe the next-to- leading order corrections, are shown in 

1 IP T'M VI H I ll (> / I = 1,7 III, , , . , , j ILLLV I'LV'J-Ill I j j 



Fig. 2a in terms of the variable r] = s/Am- — 1. The functions are split into the 



'virtual+soft' part (V+S) and the 'hard' part (H) in which the infrared log fc A (k = 1,2) 
singularities of the (V+S) contribution are absorbed so that the limit A — > can be carried 
out smoothly. The Sommerfeld rescattering contribution, due to the exchange of Coulomb 
gluons between the slowly moving gluinos in the final state, leads to a singularity ~ itot s / 
near the threshold, which compensates the phase space suppression. In addition, soft gluons 
give rise to a logarithmic enhancement near threshold, 

lT r = C' { + ^ >°g 2 (8/3 2 ) - MM 2 )} M 



with 



(0)thr = 277T P 
99 - g 4 



For large parton energies, the t- and -u-channel exchanges of gluons generate an asymptoti- 
cally constant cross section a ~ a^/m?. This is to be contrasted with the scaling behavior 
&LO ~ O-l/s in lowest order. The asymptotic values of the corresponding scaling functions 
are given by /W -> 1949/(800tt) and f$ -> -177/(160tt). 
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1.2 Quark-antiquark initial states 



In addition to the annihilation via s-channel gluon exchange, also t- and w-channel squark 
exchanges contribute to the process qq — > gg, Fig. lb. To lowest order the cross section is 
given by the function 
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where we have introduced the abbreviation m 2 _ = m? 

equal gluino and squark masses due to the negative interference between the three diagrams. 



-m~. The LO cross section is small for 



This is also the case in next-to-leading order. The scaling functions and f q q associated 
with the next-to-leading order corrections are displayed in Fig. 2b. At the threshold, the 
cross section is non-zero due to the Sommerfeld enhancement, 



r{l)thr 
Jqq 



A0)thr 
Jqq 



+ 



16/3 37T 



log 2 (8/? 2 



41 

12^2 



log(8/3 2 



(6) 



with 



Jqq — o 




m? — m?" 
m? + m? i 



Since any i- or u-channel gluon exchange diagrams are absent, the scaling functions vanish 
for asymptotic energies. 



1.3 Quark-gluon initial states 

Gluino pairs can be produced in quark-gluon collisions only in higher orders, f^J = 0, 
[cf . Fig. Id with the bottom gluon line replaced by a quark line] . The corresponding scaling 
functions and are small, in the Tevatron range, compared with the qq, gg scaling 
functions as shown in Fig. 2c. The gluon exchange in the t- and -u-channel gives rise asymp- 
totically to a non-vanishing cross section of order ct^/m? with /W — > 1949/(36007r) and 
4 1) ^-59/(240vr). 



2 Results 

The final results for the cross section pp — > ggX are presented in Fig. 3 and Fig. 4. The cross 
sections of the subprocesses have been convoluted with the parton densities in the GRV [fT5| , 
CTEQ(2pM) |14| and MRS(H) [|l5j parametrizations. The detailed analysis of these figures 
leads us to the following conclusions. 

(i) It is obvious from Fig. 3a that the theoretical predictions for the pp — > gg production 
process are greatly stabilized by taking into account the next-to-leading order QCD cor- 
rections. While the dependence on the renormalization/factorization scale Q is quite steep 
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and monotonic in leading order, the Q dependence is significantly reduced in next-to-leading 
order for reasonable variations of the scale. Close to Q ~ rrig/2 even a broad maximum 
develops. The variation of the cross section for different NLO parton parametrizations is 
very small. 

(ii) In contrast to squark-pair production, the ^-factors for gluino pair production, K = 
cnlo/clo (with all quantities in the numerator and denominator calculated consistently 
in NLO and LO, respectively), depend rather strongly on the gluino and squark masses, 
Fig. 3b, in particular in the range where squark and gluino masses nearly coincide. 

(iii) In Fig. 4 we illustrate the impact of the QCD corrections on the experimental bound 
on the gluino mass. The scale y/§ had been adopted in the experimental analyses, which 
were based on LO cross sections. We therefore compare the LO cross section, defined at this 
scale, with the NLO prediction evaluated at the scale Q = rrig. As evident from the figure, 
this increases the bound on the gluino mass by a shift between 10 and 30 GeV, depending 
in detail on m g itself and on rriq. By comparing the NLO prediction with the LO result in 
Fig. 4, it is clear that the bounds extracted by the Tevatron experiments are conservative 
and that the true bounds are likely to be higher at a level of 10 to 30 GeV. 

Acknowledgment. We thank S. Lammel for useful discussions on the Tevatron squark 
and gluino mass limits. 
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Figure 1: Generic diagrams for gluino-pair production: a) Born- level gg diagrams; b) Born- 
level qq diagrams; c) QCD and SUSY vertex corrections; d) Bremsstrahlung diagrams, and 
qg parton processes if the bottom gluon line is replaced by a quark line. 
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Figure 2: The scaling functions for gluino-pair production in a) gg, b) qq and c) qg collisisons. 
The notation follows eq. (Q) with 77 = s/4m~ — 1; mass parameters: m~ g = 250 GeV and 
nig = 200 GeV. 
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Figure 3: Total cross section for gluino-pair production at the Tevatron energy y/s = 1.8 
TeV. a) Dependence on the scale Q for the LO and NLO predictions, and sensitivity to 
different parton densities; mass parameters: m g = 250 GeV and rriq = 200 GeV. b) K 
factors for the scales Q/rrig = 2,1,1/2, corresponding to upper, middle and lower curves, 
respectively; GRV parton densities. 
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Figure 4: Total cross section for gluino-pair production at the Tevatron energy y/s = 1.8 TeV. 
Dependence of the cross section on the gluino mass; LO with EHLQ parton densities and 
Q = and NLO with GRV parton densities at the scales Q/rrig = 2, 1, 1/2 corresponding 
to lower, middle and upper curve, respectively. 
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